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Abstract 

We investigate two models for traffic flow with modified acceleration ('slow-to-start') 
rules. Even in the simplest case Vmax = 1 these rules break the 'particle-hole' symme- 
try of the model. We determine the fundamental diagram (flow-density relationship) 
using the so-called car-oriented mean-field approach (COMF) which yields the exact 
solution of the basic model with Vmax = 1- Here we find that this is no longer true for 
the models with modified acceleration rules, but the results are still in good agree- 
ment with simulations. We also compare the effects of the two different slow-to-start 
rules and discuss their relevance for real traffic. In addition, in one of these models 
we find a new phase transition to a completely jammed state. 
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1 Introduction 



The investigation of traffic flow has attracted a lot of interest from physicists in 
the last few years. In addition to more traditional methods borrowed e.g. from 
hydrodynamics, ideas from statistical physics have become more and more important. 
A comprehensive review of the different approaches can be found in |^. 

Here we use a cellular automaton (CA) P] which has been shown to reproduce at least 
qualitatively the features of real traffic 0, ^. The basic idea is to formulate a model 
in discrete space and time. The road is divided into cells which can either be empty 
or occupied by exactly one car. The cars are characterized by an internal parameter 
V = 0,1, . . . ,Vmax which corresponds to their momentary velocity. Starting from a 
given configuration at a time t the configuration at the next time step t + 1 can be 
obtained by application of four rather simple rules which are applied to all cars in 
parallel: acceleration, slowing down due to other cars, randomization and the actual 
car motion. 

The stochastic nature of the model which enters through the randomization process 
seems to be important for a realistic description of traffic. It is quite surprising that 
such a simple rule can mimic several effects which seem to be important for the 
behaviour of real traffic, e.g. overreactions at braking or fluctuations in driving. 

Due to their inherent discrete nature CA are ideal for large-scale computer simula- 
tions. On the other hand, these models can not be described easily using analytic 
approaches. In we solved the special case of v^ax = 1 exactly. We calculated the 
so-called fundamental diagram, i.e. the relation between flow (current) and the den- 
sity c of cars. The case Vmax = 1 is special since the model is particle-hole symmetric 
and therefore the fundamental diagram is also symmetric with respect to c = 1/2 
where also the flow takes its maximum value. From measurements on real traffic one 
knows, however, that the maximum flow is shifted to much lower densities. This fact 
is reproduced by the models with Vmax > 1, but unfortunately these models are no 
longer exactly solvable. 

In 10, P two models with slightly modified acceleration rules have been introduced. 
One effect of these new rules is a breaking of the particle-hole symmetry in the case 
Vmax = 1- In this paper we want to address the question whether these modified 
versions are still exactly solvable or if the exact solvability is somehow connected to 
the existence of the particle-hole symmetry. Another interesting point is a compar- 
ision of the two rules. Both are designed to mimic the same effects observed in real 
traffic, but use different ways to implement these in the rules of the CA. We therefore 
want to compare the effects of the two different slow-to-start rules and discuss their 
relevance for modelling real traffic. 



2 



2 The modified models 



In the following we will always assume Vmax = 1- The fundamental diagram of the 
basic model is symmetric (maximum flow at density Cmax = 1/2) due to a particle- 
hole symmetry . In order to get a more realistic asymmetric fundamental diagram 
even for Vmax = 1 modified models with new acceleration rules have been introduced 

Si- 

In the Benjamin- Johnson-Hui (BJH) model cars which had to brake due to the 
next car ahead will move on the next opportunity only with probability 1—ps. Each 
cell can be in one of four states. Either it is empty (state 'e') or it is occupied by a 
car with velocity v = 1 (state '1'), a car with velocity v = which is static due to the 
action of the randomization step or the slow-to-start rule (state 's'), or a car with 
velocity v = which is standing due to blocking by a car just in front of it (state 
'b'). 

1) Acceleration: All cars are assigned a velocity of '1'. 

2) Slow-to-start rule: All cars that are legitimate candidates decelerate to 's' 
with probability pg. 

3) Blockage: All cars which have no free cell in front and are thus blocked from 
moving change to state 'b'. 

4) Randomization: All cars still in state '1' decelerate to state 's' with proba- 
bility p. 

5) Car motion: All cars in state '1' move one cell ahead. 

In the Takayasu (T^) model ^ a standing car will move only if there are at least 
two free cells in front of it. More generally it will move only with probability 1 — pt 
if there is just one free cell ahead. Here every cell is one of three states, i.e. empty, 
or occupied by with velocity f = or f = 1. 

1) Acceleration: Standing cars accelerate to velocity f = 1 if there are at least 
two empty cells in front. 

2) Slow-to-start rule: Standing cars with just one free cell in front accelerate 
only with probability qt = 1 — Pt to velocity v = 1. 

3) Slowing down (due to other cars): If the distance d to the next car ahead 
is not larger than v {d < v) the speed is reduced to d — 1 [v d — 1]. 

4) Randomization: With probability p, the velocity of a vehicle (if greater than 
zero) is decreased by one [v ^ v — 1]. 

5) Car motion: Each vehicle is advanced v cells. 
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The idea behind the modified rules of the BJH and models is to mimic the delay of 
a car in restarting, i.e. due to a slow pick-up of engine or loss of the driver's attention. 
Both modifications break the particle-hole symmetry and should yield fundamental 
diagrams with a maximum at a density Cmax < 1/2. 

In order to investigate the effects of the modified rules analytically we apply the 
so-called "car-oriented mean-field theory" (COMF) This approach turned out to 
yield the exact solution for the standard rules {ps = 0, = 0) 0. It seems therefore 
worthwhile to investigate whether one also obtains the exact solution for the modified 
models without particle-hole symmetry. 



3 COMF for BJH model 

We first investigate the effects of the BJH version of the slow-to-start rule. 

We denote the probability to find at time t (exactly) n empty cells in front of a vehicle 
by Pn{t)- The density of cars which which obey the slow-to-start rule is denoted by 
Pi. As in p, Q we change the order of the update steps to 2-3-4-1. This change has 
to be taken into account when calculating the fiow f{c,p,ps)- It has the advantage 
that after step 1 there are no cars with velocity 0, i.e. all cars have velocity 1. 

The time evolution of the probabilities Pn can conveniently be expressed through 
the probability g{t) {g{t) = 1 — g{t)) that a car moves (does not move) in the next 
timestep. These probabilities are given by: 

git) = qY,Pn{t) + qsqPi(t) = q[l-P^{t)]-p,qP^{t), (1) 

n>l 

where we have used the normalization 

Y,Pn{t)+Pl{t) = l. (2) 

ra>0 

The probabilities can be related to the density c = N/ L of cars. Since each car which 
has the distance n to the next one in front of him 'occupies' n + 1 cells we have the 
following relation: 

V(n + l)P„(t) + 2Pi(t) = i. (3) 
^-^ c 

n>0 

As described in we then obtain the time evolution of the probabihties as (g = 1 —p, 

qs = 1 - Ps) 

Po(t + l) = g{t)[Po{t) + qPiit) + qsqPM (4) 

Pi(t + 1) = git)Poit) (5) 

Pi(t + 1) = g{t)\pPi{t) + {ps + qsP)Pi{t) + qP2{t)]+qgit)[Piit) + qsP,it)] (6) 

P2(t + 1) = g{t)\pP2{t)+qPs{t)]+g{t)[pP,{t) + {p, + q,p)P,{t) + qP2{t)] (7) 

P„(t + 1) = git)\pPn{t) + qPn+i{t)] + g{t)\pPn-i{t) + qPn{t)] (n>3) (8) 
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Here we are mainly interested in the stationary state {t oo) with limj^oo -Pn(^) = 
Pn- In 1^ we used generating functions to solve a similar set of equations {ps = 0). 
Here we adopt a more direct method since it is easy to see that (||) is solved by the 
Ansatz 



Pn = Nz^, {n>2) 
P9 

z = — , 
19 

where A/" is a normalization. 

Using dl-H) we can express all probabilities through Pq: 



(9) 



A = Sf^il^, (10) 

1 + PsqPo 
Pi g(l-Po) 

' = ^ = TT^' ^''^ 

Pi = ^-g.Pi, (12) 
P 

M = (i-Po-Pi-P,)i^. (13) 

The only remaining free parameter Pq can now be related to the density c by using 
d). Po(c) is given as the root in the interval (0, 1) of the cubic equation 

cplq^Pl + q\qpl{\ - 2c) +p,(l + c) + c\Pl + \qp,{\ - 3c) - 2qc + l]Po - pc = 0. (14) 
The flow is given by (with q = ^- — p) 

f{c,p,ps) = eg = c- ^ — (15) 
1 + PsqPo 

Results for the fundamental diagram and a comparison with computer simulations 
are shown in Figs. |1| and |^. Obviously, the COMF is no longer exact for ps > 0, but 
it is still an excellent approximation, especially for small ps. For fixed p (see Fig. 
|T]) the flow f{c,p,ps) is only reduced slightly for ps > 0. In addition, the density 
Cmax of maximum flow is shifted towards smaller densities with increasing ps, i.e. 
the fundamental diagram is no longer symmetric. This reflects the absence of the 
particle-hole symmetry. There is no other qualitative change of the fundamental 
diagram, in contrast to the T^ model (see next section). 



4 COMF for the model 

The probability g(t) {git) = 1 — g(t)) that a car moves (does not move) in the next 
timestep is now given by: 

g{t) = qJ2Pn{t)+qtqPiit) = q[l-Poit)]-ptqPi{t), (16) 

n>l 
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We obtain the time evolution of the probabihties as (g = 1 — p, = 1 — pt) 



Po(t + l) = gmPo{t) + qPi{t) + qtqP^{t)] (17) 
Pi(t + 1) = g{t){pt + qtp)Pi{t)+g{t)Po{t) (18) 
Pi(t + 1) = g{t)[pP^it) + qP,it)]+qgit)[P,it) + qtP^it)] (19) 

P2(t + l) = ^(t)bP2(t) + gP3(t)]+^?(t)[pPi(t) + (pt + gtp)Pi(t) + gP2(t)](20) 

P„(t + 1) = g{t)[pPn{t) + qPn+i{t)]+g{t)[pPn-i{t) + qPn{t)] (n>3) (21) 

Here we are mainly interested in the stationary state {t oo) with lim^^oo Pn{t) = 
Pn- In this case the equations (pTD - which are identical with the equations (|) since 
the modified rules do not affect distances larger than 1 - are again solved by the 
Ansatz 

P„ = Mz", (n>2) 
. = (22) 

qg 

where A/" is a normalization. Expressing all quantities through Pq and Pi we find 

Pi = -Pi = -Pi, (23) 
qg p 

g = q{l-Po)-ptqPu (24) 

= (i-Po-Pi-Pi)i^. (25) 

Pq and Pi can be related through 

PtxP^ - (1 +pgt(l - Po))Pi + (1 - Po)Po = (26) 
which yields explicitly 
1 



Pi 



1+pqtil - Po) - Vi^+PQtil - PoW - 4p*a;(l - Po)Po (27) 



2ptx 

where we have used the abbreviation x = pt + qtP = 1 — qqt- 
The relation with the density is 

y{c - 1) + (p + (1 + q)y)c{l - Pq) - (1 + p - py)cP, = (2^ 



with y = Pq + ptPi- After solving (|28D for Po(c), Pi(c) the fiow can be obtained 
using f{c,p,pt) = eg. In Figs. ^ and | we show the results of COMF and computer 
simulations. Again the COMF is no longer exact for pt > 0, but it is still a good 
approximation. For <C 1 the situation is somewhat similar to the BJH model, i.e. 
the fiow is reduced slightly and the Cmax is shifted towards smaller densities. For large 
Pt, however, one finds a qualitative change of the fundamental diagram. For ^ 1 
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the fundamental diagram has an inflection point at a density c > Cmax, in contrast 
to the case 1, where the flow vs. density relation is convex. The parameter pt 
therefore controls the curvature of the fundamental diagram. 

This can be seen most clearly for pt = 1. Here the flow f{c,p,pt) vanishes for 
c > 1/2. For c = 1/2 the stationary state consists of standing cars with exactly 
one empty site in between them (0.0.0.0. • ■ • where '0' denotes a standing car and 
'.' an empty site). All cars cannot move due to the slow-to-start rule. For larger 
densities the stationary state is essentially of the same form, but larger clusters of 
standing cars will appear. In Monte Carlo simulations starting from a mega-jam the 
relaxation into this stationary state is extremly slow, even for 'large' values of p, e.g. 
p = 0.5. Note that the COMF solution (^)-(^) predicts a phase transition to the 
completely jammed state at c = 2/3. There is, however, a second solution^ of the 
COMF equations with Pq = 1 — Pi = and P„ = for n > 1, corresponding to 
/(c,p,pt = 1) = for c> 1/2. 

The new phase transition for pt = 1 only exists for p > 0. In ^ mainly the 
deterministic case p = has been investigated. There the equilibrium state depends 
on the initial conditions and therefore the completely jammed phase for c > 1/2 could 
not been found. More detailed results for the model with pt = 1 will be presented in 
a future publication. 

5 Discussion 

In this paper we have investigated two cellular automata for traffic flow with modified 
acceleration rules. For Vmax = 1 these rules break the particle-hole symmetry of the 
model. We used the so-called Car-Oriented Mean-Field theory to obtain an analytic 
desription of these models. In contrast to the case of the basic model p the COMF 
does no longer produce an exact solution of the modified versions. 

The main effects of the slow-to-start rules are a (small) reduction of the flow and a 
shift of the location Cmax of the maximum flow towards smaller densities. Note, how- 
ever, that the reduction in the T^ model is larger than in the BJH model, especially 
for large pt- For = 1 we found a phase transition to a completely jammed phase 
with flow f{c,p,pt = 1) = near c = 1/2. This transition will be investigated in 
more detail in a future publication. 

The existence of this transition reflects the most important difference between the 
two slow-to-start rules. The BJH rule is of temporal nature (the restart probability 
depends on the number of attempts) whereas the T^ rule is of spatial nature (the 
restart probability depends on the number of empty cells in front of the car). 

In Figs. ^ and ^ we show the numbers Pi{c) and Pi{c) of cars affected by the slow- 
to-start rule in the BJH and T^ model, respectively. As expected the modifled rules 

^This solution corresponds to ( ^ ) with a '+' in front of the square-root. 
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are most effective for densities slightly larger than the density of maximum flow in 
the original model, i.e. for c > 1/2. This effect is most pronounced for small values 
of p. 

Comparing the effects of the two modifications, the model appears to be more 
interesting since the effects are larger. An interesting feature of this model is that the 
curvature of the fundamental diagram for densities c > 1/2 may vary, depending on 
the value of pt- The idealized shape of the fundamental diagram has been discussed 
controversely in the literature (see e.g. |jlO|) and there exists experimental evidence 
for fundamental diagrams with and without inflection points. The slow-to-start rule 
of T^ offers a simple explanation for these observations. 

In the present paper we only investigated the case Vmax = 1- The main modification 
introduced by the slow-to-start rules is a asymmetry between acceleration and decel- 
eration processes. Since this is already included in the basic model for Vmax > 1 we 
do not expect that the modifications might have a drastic effect for larger velocities. 
This expectation seems to be justified as preliminary results from computer simula- 
tions show. For the T^ model, however, the slow-to-start probability pt still controls 
the curvature of the fundamental diagram. 

Acknowledgment : 

Part of this work has been performed within the research program of the SFB 341 
(Koln-Aachen-Jiilich) . 

References 

[1] D.E. Wolf, M. Schreckenberg, A. Bachem (Eds.): Traffic and Granular Flow, 
World Scientific, Singapore (1996) 

[2] S. Wolfram: Theory and Applications of Cellular Automata, (World Scientific, 
Singapore, 1986) 

[3] K. Nagel, M. Schreckenberg: J. Phys. I France 2 2221 (1992) 

[4] M. Schreckenberg, A. Schadschneider, K. Nagel, N. Ito: Phys. Rev. E51, 2939 
(1995) 

[5] A. Schadschneider, M. Schreckenberg: J. Phys. A26, L679 (1993) 
[6] A. Schadschneider, M. Schreckenberg: J. Phys. A30, L69 (1997) 
[7] S.C. Benjamin, N.F. Johnson, P.M. Hui: J. Phys. A29, 3119 (1996) 
[8] M. Takayasu, H. Takayasu: Fractals 1, 860 (1993) 



8 



[9] M. Fukui, Y. Ishibashi: J. Phys. Soc. Jpn. 66, 385 (1997) 
[10] F.L. Hall, B.L. Allen, M.A. Gunter: Transpn. Res.-A 20A, 197 (1986) 



9 



p=0.5, p, =0 p=0.5, p, =0.1 




Figure 1: Fundamental diagram for the BJH model with p = 0.5 and different Ps- The 
full line is the COMF result. For comparison the results from computer simulations 
(•) are also shown. 
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p = 0.001, =0.5 p = 0.1, p, =0.5 

-] 1 1 1 1 1 1 1 1 I 1 1 1 1 1 1 1 1 r 
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Figure 2: Same as Fig. |l|, but for ps = 0.5 and different p. 
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p=0.5. pt =0 p=0.5, pt=0.1 




Figure 3: Fundamental diagram for the model with p = 0.5 and different ps- The 
full line is the COMF result. For comparison the results from computer simulations 
(•) are also shown. 
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p = 0.001, pt =0.5 p=0.1, pt=0.5 
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Figure 4: Same as Fig. y, but for ps = 0.5 and different p. 
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BJH model 




.2 .4 .6 .8 1.0 



C 

Figure 5: Pi for the BJH model for a) p = 0.001, Ps = 0.5, h) p = 0.1, Ps = 0.5, c) 
p = 0.5, Ps = 0.1, d) p = 0.5, = 1, e) p = 0.75, p^ = 0.5. 
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Figure 6: Pi for the model for a) p = 0.5, = 1, b) p = 0.5, Ps = 0.75, c) 
p = 0.75, Ps = 0.5, d) p = 0.1, Ps = 0.5, e) p = 0.001, = 0.5. 
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